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We study a five-dimensional perfect fluid coupled with Kaluza-Klein (KK) gravity. By dimensional 
reduction, a modified form of Maxwell's equation is obtained, which is relevant to the equation of 
state of the source. Since the relativistic magnetohydrodynamics (MHD) and the 3-dimensional 
formulation are widely used to study space matter, we derive the modified Maxwell's equations 
and relativistic MHD in 3-1-1 form. We then take an ideal Fermi gas as an example to study 
the modified effect, which can be visible under high density or high energy condition, while the 
\ traditional Maxwell's equation can be regarded as a result in the low density and low temperature 

. limit. We also indicate the possibility to test the state-relevant effect of KK theory in a telluric 

C^) ' laboratory. 
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I. INTRODUCTION 



O^' A unified formulation of Einstein's theory of gravitation and Maxwell's theory of electromagnetism in four- 
dimensional (4D) spacetime was first proposed by Kaluza and Klein using a five-dimensional (5D) geometry [l|, 0]- 
A free test particle in 5D KK spacetime shows its electricity in the reduced 4D spacetime when it moves along the 
fifth dimension. Moreover a 5D dust field coupled with KK gravity can curve the 5D spacetime in such a way that it 
provides exactly the source of the electromagnetic field in the 4D spacetime after the reduction . In this paper we 
study the coupling of a 5D perfect fluid with KK gravity. It turns out that the 5D Einstein's equation with a source 
gives a modification of Maxwell's equation (see in section II) , which can show its state- relevant effect on high-density 
. or high-temperature condition (see in section IV). Thus this effect provides intriguing possibilities for the experimental 

I \ test of the KK theory. Note that the KK theory which we are considering is purely classical. The KK theory has also 

I been studied from the particle physics point of view (see e.g. 0|), which can show its effects in high energy scale of 
TeV . Whereas the energy scale needed for testing our state-relevant effect is only around keV . 

In order to reveal the physical implication of the modification more clearly, both the modified Maxwell's equations 
. !^ and the corresponding general relativistic MHD are reformulated in 3+1 form (see in section III). The formalism is also 
' useful for evolving numerically a relativistic MHD fiuid in a spacetime characterized by a strong gravitational field. 
^ . Taking an ideal Fermi gas as an example, the modification term is studied as a function of degeneracy and temperature 
parameters (see in section IV). Moreover the modification terms for different components of the perfect fiuid may 
be different. This would result in a net charge excess in high-temperature plasma. Recall that the the electrical 
neutrality of atoms and of bulk matter has been examined precisely by a number of experiments However, in 
those experiments, either the objects considered are not ionized, or the ions in the objects cannot be regarded as a 
perfect fluid. Therefore, these experiments cannot provide definite opposite evidence to the classical KK theory since 
they do not satisfy our premise. But they do cast some doubts on it. Taking account of the state-relevant character, 
we suggest high-temperature plasma in earth laboratory or dense-matter white dwarf in outer space as candidates to 
test the possible effects of the modifled Maxwell's equation. 
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II. KALUZA-KLEIN GRAVITY COUPLED WITH 5D PERFECT FLUID 



Using a 5D geometry, Kaluza and Klein proposed a unified formulation of gravity and electromagnetism in 4D 
spacetime P, The original KK theory assumed the so-called "CV//OT(fer Condition" , which means that there exists 
a space-like killing vector field on the 5D spacetime {M,gab) 0-01 • Note that the abstract index notation 9] is 
employed throughout the paper and the signature of the five-metric is of the convention (— , +, +). In addition, 
Kaluza also demanded that is normalized, i.e., 

<P = gabCt = 1 ■ (1) 

Later research shows that the Ansatz ([1]) may be dropped out and the 4> may play a key role in the study of 
cosmology [10|, [IJ, Being an extra dimension, the orbits of are geometrically circles. The physical 

consideration that any displacement in the usual "physical" 4D spacetime (denoted as M) should be orthogonal to 
the extra dimension implies that the "physical" 4D metric should be defined as 

9ab = gab - <f>^^Lib , (2) 

and the projection operator onto M is 

9% - 9% - 0"'r6 • (3) 

For practical calculation, it is convenient to take a coordinate system {z*^ — (a;^,?/)|^ = 0,1,2,3} with coordinate 
basis (gm)" = {(e^)", (65)°} on M adapted to i.e., (65)" = (^)" — C"- Then the 5-metric components qmn take 
the form 

= ( ) , (4) 

where g^^ = (j)B^. So, locally, the "physical" spacetime can be understood as a 4- manifold M with the coordinates 
{x'^} endowed with the metric gab- The whole theory is governed by the 5D Einstein-Hilbert action 

Sg^~^ I d^xdy^gk . (5) 

Suppose the range of the fifth coordinate to be < y < L and denote k — k/L. Let 5^ — JA^^ p — 2k, then equation 
([5]) becomes a coupling action on M as 



Sg = 



jj^x^g^(^-^R+\<t,Fab{A)F'^\A)^ , (6) 



where R is the curvature scalar of gab on M and Fab{A) = 2d\^aAb]- Thus, it results in a 4D gravity gab coupled to 
an electromagnetic field Aa and a scalar field (f). It is clear that, under the Ansatz ([T]), 5D KK theory unifies the 
Einstein's gravity and the source-free Maxwell's field in the standard formulism. 
Now we consider a 5D perfect fluid 

Tab (P + i^WaVb + Pgab ■ (7) 

The 5-velocity V'^ can be projected onto the "physical" spacetime (Af , gab) as 

= g^V"^ = V^^ie^r - {B^Vn{e5r ■ (8) 
Note that we have V"'Va = — 1, hence it is easy to show that 

V^^V^ = gabu'^v!' = gabV^'V" = -1 - ^ , (9) 
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where Q = V5 = .gsaV'" represents the electric charge in M The energy- momentum tensor can be projected on M 
as Tab = ga9b^Tcd- In order to obtain the observed 4D energy-momentum tensor Tab on M, we have to integrate Tab 
along the extra dimension. In the light of ([5]) and (O we obtain 

Tab = (m + P)VaVb + Pgab , (10) 

where 



P 



= P\/<i>L, 



Ua 



(11) 



It is clear that Tab is the energy-momentum tensor of a 4D perfect fluid in M, where /i and p are respectively the 
energy density and pressure density observed by a comoving observer in M . 
We now consider the reduction of 5D Einstein's equation 



which is equivalent to 



Rab ~ -^cjabR = kfab , (12) 



Rab = HTab - ^f^gab) ■ (13) 



It is not difficult to show from Eq. ([5]) that the components of the 5D Ricci tensor Rab can be expressed as 
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i?55 = ^kcj^^F'^PP^p _ i + ^ (V^^) , (14) 



= ((^V-F^. + ^F^.V^ 



-B„ 



(15) 



R^. = R^.. - k<t>F^^,F,, - ^V^V,0 + -L(v^</,)V,0 

+B^B, (^kcj^'^F'^PP^p - iv"V<,</.-K ^(V"</))V,0^ 
+ |b^((/)V'^F,, + ^F.^V"^) 

+ ^B,i4'V^F^, + ^F^aV''^) , (16) 

where Va is the 4D covariant derivative operator associated with gab- Substituting Eq. (|14p into Eq. (jl3|) . we obtain 
a coupling equation for the matter fields as 

h^^F^^Fab = y^V'^V.v/^+ fcV^. (l - + ^V^P^^^ . (17) 

Substituting Eq. (fT51) into Eq. (|13p and using Eq. PTI) . we obtain an electromagnetic field equation with source as 

'/'V'F.fc + ^F^fcVV = 7(1 + -)Ja , (18) 
2 /Lt 
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here we have defined 7 = \/(TTW)/W+W) , P = ^^{^^ga) J" = P'"" 3- Substituting Eq. ^ into Eq. JTSl) 
and using Eq. (fT7|) and (fT8|) . we obtain a 4D Einstein's equation with source as 



Gab = ^ (('^ + P)^'^^'' + 3abP + <l>^^HF^'Fbc - ^9abF^''F,d) 



(19) 



where Gab is the Einstein tensor of gab- More generally, if the 5D perfect fluid consists of m components, Tab then 
reads 



Tab^Y.{(Pn + f'vWa{v)Vb{v)+P,i9ab) ■ (20) 

ri=l 

By similar calculations, Eqs. P?)) - p^ become respectively 

h^^F-'Fab = V^V'^Vav^+A:v/^f^A*('7) fl ^'^ 



r,= l 



3(Q(r7)2 + </.) 



'/'V^f^afc + ^Fab^'c^ = pliv) (l + ^) Jaiv) 



= 71 (^i^^^^'^^ +Piv))va{v)vbiv) + 9abP{v)) + <t>^'\F^-Fbc - -gabF^'^F.a 



(22) 



-^(ffafcV^V.V^- VaV^y^)") . (23) 



It is interesting to see the results when 0=1. Eqs. (fT7|) -(fT9 |) become respectively 

wrm) + ''m^) ' ^''^ 

^'Fab = (1 + -) Ja , (25) 

G.,=fc(Tir'^)+rir^), (26) 

where T'^^""^'' = (a* +p)vaffc + pgab and = F^^Fbc — jgabF^'^Fcd are respectively the usual energy-momentum 

tensors of 4D perfect fluid and electromagnetic field. Eq. (PSI) is the standard 4D Einstein's equation while Eq. P5|) 
is not the same as the standard 4D Maxwell's equation V^Fab = Ja- The new term (1 + ^) brings an effective charge 
which can be considered as a state-relevant effect, as will be discussed later. We thus call Eq. ((25l) the state-relevant 
Maxwell's equation. 



III. MAXWELL'S EQUATIONS AND RELATIVISTIC MHD IN 3+1 FORM 



Since relativistic MHD is widely used to study space matter and 3D formulation is frequently applied to dealing 
with specific issues [14| . we now derive the modified results of Maxwell's equations and relativistic MHD in 3-1-1 form. 
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The 4D spacetime M is foliated into a family of non-intersecting spacelike three-surfaces E, which arise, at least 
locally, as level surfaces of a scalar time function t. The spatial metric ^ab on the three-dimensional hypersurfaces S 
is induced by the spacetime metric Qab according to 

lab = gab + naUb , (27) 

where n° is the unit normal vector to the slices and thus Ua — — aVat- Here the normalization factor a is called the 
lapse function. The time vector i° is dual to the foliation 1-form V at and can be decomposed as 

= an" + (3" , (28) 

where the shift vector /3 is spatial, i.e., rta/?" = 0. Since the extrinsic curvature Kab of S can be written as 

Kab = -VaHb - Uattb , (29) 

where aa = n^WbTia, the divergence of satisfies 

V,n° = -K , (30) 

here K is the trace of Kab. 

Firstly we write the modified Maxwell's equations in 3+1 form. The Faraday tensor can be decomposed as 

pab ^ ^^apb _ ^bpa ^ ^abc^^ ^ (3^^ 

where E°' and are the electric and magnetic fields observed by a normal observer n". Both fields are purely spatial, 
whereby 

E^Ua = and B^Ua = , (32) 
and the three-dimensional Levi-Civita symbol eabc is defined by 

_ abc _ dabc /oo\ 

eabc = n edabc or e = n^e . (33) 

The electromagnetic current four-vector J" is decomposed as 

J" = n>e+r, (34) 

where Pe and are the charge density and 3-current as observed by a normal observer Note that j" is purely 
spatial, i.e., j"'na — 0. With these definitions, the modified Maxwell's equations (|17p . p8)) and W[aFbc] — can be 
cast into 3-1-1 form as 

k<j>\B^-E^) = ^{D-Da^^{a-\dt- Cfi) f^ 
+Ka-\DaV^){D^ liia)idt - Cfj)^) 

DaE" = 0-1(7(1 + ^)p,-^£;-i?,0), (36) 
/i 2 

(j)£tE'' = (l){aKE'' +CpE'' + Db{aBc)) 

-a7(l + + ^a{e-'-BbD,q^ - E-a-\dt - Cp^ , (37) 

DaB'' = , (38) 
CtB" = -e''^''Db{aEr) + aKE'' + CpB" , (39) 
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here denotes the Lie-derivative along /3° and Da is the covariant derivative operator associated to ^ab- Note that the 
Lie-derivative of a spacehke tensor A'^^'^^ along is defined conventionally as £.sA°'c--d = 7e ' ' ' 7/ 7c ' ' ' Id^sA'^gJj,, 
and we write as Cs for short. Note also that the formula n^'VijUa = = Dalna is used in the above calculation [15j . 



If one considered a 5D perfect fluid consisting of m components, the terms 7(1 + j^)pe and 7(1 -|- in Eqs. ((36l) 

and ((37| would be replaced by + j^)Peiv) and X^rJLi 7('/)(l + When 0=1, one can see from 

Eq. ([55]) that pe = (1 + p/ fi)pe is the effective charge density serving as the source of the electric field. This effective 
charge density is state-relevant, i.e., it is dependent on p/ p. Its significance will be discussed later. 

Secondly we rewrite modified relativistic MHD in 3+1 form. Note that the total energy- momentum tensor in M 
can be read off from the right hand side of Eq. (|19p as 

T'-'^T^L^^+f^^^^^+T^^) , (40) 

where 

f^,l^ ^ ^^'^T^L) = <ly"^{F^'F\ - ig'^'F^'^F.rf) , (41) 
Tf^) ^ -i(5"''V^Vc\/^-V'^V^y^) . (42) 
It is straightforward to see that [l3| 

VfeTfl) = r^-^^'Fbc . (43) 
In the light of Eqs. (I17I)-(I321), we obtain the 3+1 form 

-E^'DbE'' + e°''"'B^a-^{Ct - Cf3 - aK)Eb 

+a-^BbD''{aB'')^a-^BbD\aB'') , (44) 



2 (B^ + i?^)(i?> + n"«-i(at-/:0)0) 



-{E'^E'' + B''B'')D, 



+EMn-e''-^Dhcj) + e'^-'*a-i(at - £^3)0)) , (45) 

here E^ = EaE"" and B^ = BaB"' . Recall that the relation between three-dimensional Riemann tensor '^R^bc 4D 
one R'^^^^ reads 

'Rtc = lllililtKfi - '^K^laK^. (46) 

Using 

(V,V, - V,V,)V^v^ = -i?,JV'^y^, (47) 



a lengthy but straightforward calculation gives 

1 

-(i^V)(«"'(A - Cp)K% + ^"'-^fe 
-D'^Dh In a + [D" In a)Dh In a) 



+n"(a-i((9t - Cp)^){~K^'Kbc + i^'^'i^'b Ina + 
+a-\dt - /:/3)X) + {DbK^')D,(l> - {D''K)Db(j) 

-KbciD''^)D'lna)), (48) 



where "^R is the three-dimensional Ricci tensor and a = a^Ua- For a perfect fluid, the energy-momentum tensor 



T^J" can also be written as 



T^l,,)=phv-v''+pg'^'' , (49) 



where p is the rest-mass density as observed by an observer co- moving with the fluid v"^, p is the pressure and h the 
specific enthalpy 

h = l + e+p/p. (50) 
Hence one has p, = {1 + e)p. The local conservation of the 4D Einstein tensor G"'' leads to 

Vf,(T'^Vv^) = . (51) 
We assume that T^^-^ does not contribute to the number of baryons. Thus we have the conservation of baryons as 

Va(p«'^) = 0, (52) 

which is decomposed into 3+1 form as 

Daipv") + a-\dt ~ li(i){pW) ~ pWK -f pi^Dalna = 0, (53) 

where -D" = — Wn"^. The equation for the conservation of energy is obtained by contracting Eq. ([51]) with ri^ as 

H = a-\dt~ Cp)p-Wpiv^Dah + Wa-\dt- Cf3)h) 
-ph{Wa~\dt - Cp)W + v^DaW 

+Wi''Da\na~ Kabi'^v''), (54) 

and the Euler equation is obtained by projecting Eq. (j5ip onto S as 

phv^Dbv'' = ph{2WK''''ib - W^D" In a) - phWa-\£t - Cp)ij'' 

-pi^iWa-^dt - Ci3)h + i^Dbh) - D^p - M", (55) 



where 



+ ^(^2 ^ E^)a-'idt - Cp)cp + ^e'^'^'EaBbD^c^ 

+ ^^(a-'(5t - Cp)c^{-K''''Kab + D'^Da Ina + a^ 
+a-\dt - C[3)K) -f {DaK''')Dbcb - {D-K)Da(b) 
-^{p + p)v'Db^ + -^{DaCp){a-\Ct - Cp)D-^ 

-KabD''^) + l.a-\-W^p+p)+p^UD''Day/^ 
Zcp k 

+Ka-\dt - Cp)^+{Da^)D''\na)){dt - Cp)4>, (56) 



and 
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+e-'''EtB,a-\dt - Cp)ci,) - (jy^'^E'^D^E" 

+(j>^/^a-\e'''"'Bc{Ct -Cp- aK)Eb + BbD''{aB'') - BbD^iaB")) 

+ f + a-\DbK^' - D^K){dt ~ Cp)4>) 

- {D^(t)) {a-^{Ct ~ Cp)K\ + K'^'^Kbc + D'^Db In a + {D" In a)Db In a) 

+Ka-\dt - Ci3)y/^+{Dby/^)D' In a) 

-{a-^dt - Cp)4,)D-{a-\dt - Cf3)V^)) , (57) 

here {a-^{dt - Cf3))^^/^ denotes a'^idt - Cp){a-^{dt - Cp)^/4>). Note that Eqs. ^ and ^ comprise the 

basic formulas for the modified relativistic MHD in three-dimensional form. In the special case ^ = 1, we have 



H = -KE' + a-'EaiCt - Ci3)E'' ^ a-'e^'^EaDbiaBc), (58) 
M" = -E''DbE'' + a-^e'''"'Bc{Ct-Cp-aK)Eb 

+a-^ {BbD''{aB'')- BbD\aB'')) , (59) 

which accord with the conventional form uA 



IV. DISCUSSION ON THE EFFECTIVE CHARGE 

From the state-relevant Maxwell's equation (|25p and its 3+1 form Eq. ([55]) . we can see that pe = (1 +p/ii)pe is 
the effective charge density of a perfect fluid. Such effective charge density is relevant to the equation of state of the 
fluid and its rationality should be carefully checked. In this section, we adopt an ideal Fermi gas as an example to 
see under what condition can the modified term p/ ji show visible effect. We employ in terms of the dimensionless 
degeneracy and temperature parameters 

where jl is the chemical potential, m is the mass of the fermion, and kB is the Boltzmann constant. The gas is 
degenerate for 77 ^ while nondegenerate for 77 ^ 0. On the other hand, the gas is extremely relativistic for /? 1 
while nonrelativistic for /3 <C 1 [itI . The zero of energy for the particles is chosen so that the thermodynamic potential 
reads 



n = -VkBT I ^In 



/l3 

where p is the momentum, g is the statistical weight, and 



cxp 



kuT 



(61) 



e = [pc)"^ + (mc^)^ — mc^ 



(62) 
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FIG. 1; Modified term as a function of degeneracy parameter ri and relativistic parameter /3. 



is the kinetic energy. The number density n, pressure p, and internal energy density E (per volume) of an ideal Fermi 
gas are respectively 

n = if/33/2 [F,/2iv,l3)+PF3/2{v,P)] , (63) 



p — mc K(3 



(64) 



E = TOc^K/?^/^ [i^3/2(r/,/3) +/3F5/2(?7,/3)] , (65) 
where K — 4:^/2ng{mc/h)^, and the Fermi integral is 

Fkiv.P)^ ^ (fc>-l). (66) 



Then we have 

P 



E + nmc^ 3 + 3Fi/2(?7,/3)/(2/3i^3/2(r7,/3)+/32F5/2(r?,/3))- 



(67) 



The relation between p/fj, and {rj^P) is shown in Fig. [T] One can see clearly that both degenerate and relativistic 
conditions can lead to the value of p/fi comparable to 1/3 (which is the value oi p/n for radiation). 

Now we study the two kinds of conditions respectively. For a nondegenerate ideal Fermi gas (for example rj = —30), 
the value of p/fi is drawn from nonrelativistic {/3 — 0) to relativistic {/3 — 2) regime in Fig. O It is obvious that 
we need not to go to extremely relativistic condition since p/fi is already close to 1/3 when (3 = 2. In the specific 
calculation for an electron gas, we set p//i = 0.002 when fc^T = 1 keV. This result indicates the possibility to test 
the theory in earth laboratory. For a non-degenerate electron gas at T = 273A', one may estimate the modification 
term as p//i ^ fc^T/meC^ ^ 10^^. On the other hand, in the experiments on the equality of the electric charges 
of proton and electron, these charges in a conductor are found to be equal within 10~^^ or better (see e.g. jl8|). 
However, the proton system in a conductor cannot be seen as a perfect fluid and hence does not satisfy our premise. 
Hence the effective charges of protons in a conductor cannot be directly obtained by our modified equations. So, 
those experiments are not in severe contradiction with the KK theory. For similar reason, the experiments reported 
in Ref.l5| cannot provide definite opponent evidence to the KK theory either. But this kind of experiments do cast 
some doubts on the classical KK theory. Note that both the electron system and the ion system could be regarded as 
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FIG. 2: Modified term as a function of relativistic parameter (3 in nondegenerate condition rj = —30. 
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FIG. 3: Modified term p//i as a function of degeneracy parameter -q in nonrelativistic condition /3 = 10 



perfect fluid in high-temperature plasma. In a thermal equilibrium state the electron and ion in a plasma have the 
same temperature. Hence they would have different values of p/fi. Actually the value of p/fi for ion is much smaller 
than the one for electron when ksT takes value from keV to MeV. It turns out that the two important physical 
parameters for the description of plasma-Debye length and plasma frequency [19] have to be modified in our 5D 
theory as 



Xd - 



and 



Uee'^il+Pe/He) 



rrieeo 



(68) 



(69) 



where Ue is the number density of electron and eo the permittivity of vacuum. Since the electromagnetic wave whose 
frequency is lower than Wj, will be reflected while others can transmit through the plasma, the plasma frequency can 
be measured accurately |20| . Therefore it is possible to test the prediction from the 5D KK theory in earth laboratory. 
For a degenerate idea Fermi gas, the relation between and j] is demonstrated in Fig. [3l Recall that the white 
dwarf is known to resist the gravity by an electronic degenerate pressure. It is also possible to test the 5D theory by 
certain relevant phenomena in outerspace. 
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Note that the vacuum polarization in quantum electrodynamics (QED) also leads to an effective charge of a point- 



like particle 2l|, . So the effective charge viewpoint does not merely come from the KK theory. For the Fermi gas in 



KK theory, the larger the density and the temperature, the larger the effective charge factor Pe/ Pe, which approaches 
to 4/3 as a limit. Whereas for QED, the higher the energy scale (or shorter distance), the larger the effective charge 
Ccff/e, which approaches to infinity as a limit. Therefore the state- relevant Maxwell's equation and QED give similar 
results of larger effective charges. However, the state-relevant effect in KK theory is a pure classical effect due to the 
extra dimension of spacetime, whereas the QED effect is a quantum effect irrespective of any extra dimension. So one 
does not expect them to be the same. It is easy to distinguish the two effects by comparing their characters. 

In Summary, the coupling of 5D perfect fluid to KK gravity is fully studied. The 4D effective equations of this 
5D coupling system are derived. In particular, the modified Maxwell's equation which is relevant to the equation 
of state of the source is obtained. To facilitate applications, we also derive the 3-1-1 form of the modified Maxwell's 
equations and the relativistic MHD. It turns out that the effective charge density in the KK theory can be written 
as Pe = (1 + pI p)pe- Moreover, using an ideal Fermi gas model, we study the modification term p/ p as a function of 
degeneracy parameter and the relativity parameter (3. It reveals that the traditional Maxwell's equation is the low 
density and low temperature limit of the state-relevant Maxwell's equation. We thus indicate the possibility to test 
the state-relevant effect both in earth laboratory and in astrophysical phenomena. 



Acknowledgments 

We acknowledge the valuable discussions with Lingzhen Guo, Wenan Guo, Zhi-Qiang Guo, Bin Wu and Ren-Xin Xu. 
This work is supported in part by Hui-Chun Chin and Tsung-Dao Lee Chinese Undergraduate Research Endowment 
(Chun-Tsung Endowment) at Peking University, by NSFC (Nos. 10675019, 10421503, 10575003), and by the Key 
Grant Project of Chinese Ministry of Education (No. 305001). 



[1] T. Kaluza, Sitzungsber. Preuss. Akad. Wiss. Phys. Mat. Klasse 1921, 966 (1921). 
[2] O. Klein, Z. Phys. 37, 895 (1926). 

[3] Y. Ma and J. Wu, Int. J. Mod. Phys. A19, 5043 (2004). 

[4] Y. A. Kubyshin, Arxiv preprint; h ep-ph/01 11027 

[5] C. S. Unnikrishnan and G. T. Gillies, Metrologia 41, S125 (2004). 

[6] M.J. Duff, B.E.W. Nilsson, and C.N. Pope, Phys. Rep. 130, 1 (1986). 

[7] J.M. Overduin and P.S. Wesson, Phys. Rep. 283, 303 (1997). 

[8] X. Yang, Y. Ma, J. Shao, and W. Zhou, Phys. Rev. D68, 024006 (2003). 

[9] R.M. Wald, General Relativity (The University of Chicago Press, 1984). 
[10] J-P. Uzan, Rev. Mod. Phys. 75, 403 (2003). 
[11] P.J.E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003). 
[12] I.K. Wehus and F. Ravndal, Int. J. Mod. Phys. A19, 4671 (2004). 
[13] N. Mohammedi, Phys. Rev. D65, 104018 (2002). 
[14] T.W. Baumgarte and S.L. Shapiro, Astrophys. J. 585, 921 (2003). 
[15] K.S. Thorne and D.A. MacDonale, Month. Not. R. Astr. Soc. 198, 345 (1982). 
[16] C.W. Misner, K.S. Thorne and J.A. Wheeler, Gravitation, (Freeman, New York, 1973). 
[17] J.A. Miralles and K.A. Van Riper, Astrophys. J. Suppl. S 105, 407 (1996). 
[18] R.W. Stover, T.I. Moran and J.W. Trischka, Phys. Rev. 164, 1599 (1967). 
[19] J.A. Bittencourt, Fundamentals of Plasma Physics, (Springer, 2004). 
[20] R. Xu, Introduction to Astrophysics, (Peking University Press, 2006, in Chinese). 
[21] W. Greiner and J. Reinhardt, Quantum Electrodynamics, Second Edition, (Springer, 1994). 

[22] M.E. Peskin and D.V. Schroeder, An Introduction to Quantum Field Theory, (Addison- Wesley Publishing Company, 1997). 



